Polar surface engineering in ultra-thin MgO(lll)/Ag(lll) — possibility of 
metal-insulator transition and magnetism 
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A recent report [Kiguchi et al, Phys. Rev. B 68, 115402 (2003)] that the (111) surface of 5 
MgO layers grown epitaxially on Ag(lll) becomes metalhc to reduce the electric dipole moment 
raises a question of what will happen when we have fewer MgO layers. Here we have revealed, first 
experimentally with electron energy- loss spectroscopy, that MgO (111) remains metallic even when 
one-layer thick, and theoretically with the density functional theory that the metallization should 
depend on the nature of the substrate. We further show, with a spin-density functional calculation, 
that a ferromagnetic instability may be expected for thicker films. 

PACS numbers: 79.60.jv, 73.20.-r, 68.55.-a 



The problem of polar surfaces, i.e., surfaces of ionic 
crystals where all the outermost atoms are either anions 
or cations, is of fundamental as well as technological in- 
terest. A textbook description is that a polar surface, 
which would accompany a spontaneous appearance of a 
macroscopic electric field, should be unstable, and the 
surface usually reconstructs itself to wash out the elec- 
tric field. While polar surfaces have been envisaged to 
be a unique playground for catalysis, unusual adhesion, 
etc., due to their high reactivity, one distinct feature is 
the possibility of metallization of the surface layer that 
would imply a purely two-dimensional metal, on which 
we focus in the present study. 

The physics behind the metallization of polar surfaces 
is as follows. A slab having polar surfaces is equivalent 
to a stack of dipole layers, which produce the bulk elec- 
tric field. Since this state has an infinite surface energy, 
a sheet of compensating charges must be formed at the 
top and bottom surfaces. The system can accomplish 
this in two ways: either the surface is reconstructed to 
introduce, e.g., systematic vacancies, or a charge trans- 
fer may occur across the bulk and the surface without 
reconstructions. In the latter case the sheets of surplus 
charges imply that each surface becomes metallized. 

The proposal of the metallized polar surface has long 
been investigated theoretically^. Specifically, the unre- 
constructed MgO (111) surface, which consists of a top 
0^~ surface and a bottom Mg^+ surface (type-HI sur- 
face in the classification by Tasker^), a quantum me- 
chanical calculation has been performed with the dis- 
crete variational(DV)-XcK method^ for a finite cluster^, 
a semi-empirical Hartree-Fock^, or a first-principles den- 
sity functional calculation^'^ for a slab. The metallization 
has been confirmed in these studies. 

On the other hand, the problem has a long history 
of experimental investigation as well, but only recently 
do we have definite results. Namely, various attempts to 
grow polar surfaces of rocksalt-structure compounds with 



unreconstructed polar surfaces have proved to be a dif- 
ficult task. It has in fact been shown theoretically that 
reconstructions cannot be avoided unless we introduce 
hydroxylation^'^ or adsorption of metals^'^'^. If we turn 
to oxide surfaces, they are even more difficult to prepare 
than ionic crystal surfaces, because oxygen deficiencies 
tend to occur and hinder spectroscopy for oxide surfaces 
even when non-polar. Now, recent advances in fabrica- 
tion techniques for ultrathin oxides^^'^^ is enabling us to 
prepare atomically-controlled oxide single crystal films, 
but the unreconstructed polar surfaces remained a chal- 
lenge. 

Recently, however, three of the present authors^^ have 
succeeded in growing an MgO (111) 1x1 surface by al- 
ternate adsorption of Mg and O2 on Ag(lll), which 
has enabled them to have a unreconstructed polar sur- 
face (which turned out to be metastable, transforming to 
more stable ones after annealing). The electronic struc- 
ture of 5 MgO layers (or 10 monolayers in another nomen- 
clature) of MgO (111), as probed with electron energy- 
loss spectroscopy (EELS) and ultraviolet photoemission 
spectroscopy (UPS), shows that the surface is indeed un- 
reconstructed with a nonzero density of states at the 
Fermi energy {Ep) indicative of a metal. 

This raises an interesting question of what should hap- 
pen to the metallization when we have only few mono- 
layers: As the film thickness decreases the macroscopic 
electrostatic potential due to the surface charge should 
decrease, so the competition between the metallic state 
and the insulating state accompanied by an electric field 
should become more subtle and interesting. The present 
study exactly addresses this problem. 

Here we show, first experimentally, that even one layer 
of MgO (111) on Ag(lll), as probed with EELS, is metal- 
lic. We then theoretically show with an ab-initio calcu- 
lation that the system is indeed metallic. We further 
predict that the nature of the surface should depend on 
the nature of the substrate, i.e., the surface should be 
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insulating for substrate with larger lattice constants. We 
finally show, with a spin-density functional calculation, 
that the surface, whose local density of states is large at 
the surface for thicker films, can exhibit a ferromagnetic 
instability. 

Let us start with the experimental result. As for the 
sample preparation we follow Ref.^^. Namely, the MgO 
film was grown by alternate adsorptions of Mg and O2 on 
the Ag substrate with temperature of 300K, where the 
whole experiments were performed in a ultrahigh vac- 
uum. The formation of a clean 1x1 MgO film on the 
substrate is confirmed from a sharp reflection high energy 
electron diffraction (RHEED) pattern, while the Auger- 
electron spectroscopy detects no contamination. 

The electronic structure is then probed in situ with 
EELS. The result is displayed in Fig.l for one layer of 
MgO(lll) on Ag(lll) as compared with MgO(lOO)^^ of 
the same thickness grown on Ag(lOO). We can immedi- 
ately see that the MgO (111) has a long and substantial 
tail for the energy loss ^ 2.5 eV in a marked contrast 
with the result for MgO (100) (while the peak around 4 
eV, visible in both cases, originates from the surface plas- 
mon of Ag substrate). The qualitative feature observed 
here for one layer of MgO on Ag is quite similar^^ to those 
of five layers of MgO on Ag, reported in Ref.^^. So we 
conclude that one layer of MgO (111) on Ag already has, 
surprisingly enough, a nonzero density of states around 
Ef suggestive of a metallic surface. 
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FIG. 1. EELS result, with the primary-electron energy of 
60 eV, for the one MgO(lll) layer on Ag(lll) as compared 
with that for MgO (100) on Ag(lOO). 



We have then performed a first-principles electronic 
structure calculation in the framework of the general- 
ized gradient approximation (GGA) based on the density 
functional theory. We adopt the exchange-correlation 
functional introduced by Perdew et al.^^ We employ 
ultra-soft pseudopotentials^^'^'' for Mg and O, and a 
norm-conserving soft pseudopotential for Ag, both in sep- 
arable forms. The cut-off energy of the plane- wave ex- 
pansion for the wave function is taken to be 25.0 Ry. The 
atomic configurations and the corresponding electronic 
ground states are obtained with the conjugate gradient 
scheme^^. 

The optimized lattice constants in the bulk obtained in 
the present calculation are 4.11 (against the experimen- 



tal 4.09) A for Ag and 4.22 (4.21) A for MgO. We then 
introduce a slab model, where we put a (lll)-directed 5 
Ag layers sandwiched from top and bottom by MgO lay- 
ers. The structural optimization is done as follows. We 
put Mg atoms on the hollow sites of the outermost Ag 
atoms (see Fig. 2), since the total energy is found to be 
lower by ~ O.leV per unit cell than when we put them 
on the atop sites. The bulk lattice constants of Ag and 
MgO are slightly different (4.09 against 4.21), and we as- 
sume in the calculation that the lattice constant of the 
Ag substrate does not change at the interface. So we 
take the same unit-cell size as the bulk of Ag in the xy- 
plane, and optimize the positions of the atoms. The size 
of the whole slab in z direction is set to be large enough 
(47.46 A). Then we put oxygen atoms to have one layer 
of MgO (111). We repeat this for > 2 layers for multi- 
layer cases. While larger unit-cell sizes in xy would be 
necessary to describe reconstructed surfaces, we take 1 
X 1 unit cell here since we focus on the unreconstructed 
case. The reasoning is as follows. The unreconstructed 
surface is experimentally observed as the 1x1 RHEED 
pattern, but the state is metastable. So if we want to fo- 
cus on such a state we can perform a density functional 
calculation in the prefixed unit cell. 




FIG. 2. Top (a) and side(b) views of the atomic configura- 
tion considered in the present study. The unit cell is indicated 
by white lines. 



Let us start with the result for three MgO layers on 
Ag. Figure 3 shows the band structure along with the 
local density of states(LDOS). The LDOS at Ep^ plotted 
in Fig. 3(c), is calculated by \(j)i{x^y, z)\'^ , where the 
summation is taken over the eigenstates (labelled by i) 
having the eigenenergies ^ — 0.125 < Ei < £^ + 0.125 eV. 
The LDOS does not change significantly when the en- 
ergy window (0.25 eV here) is changed to 0.5 eV. The 
number of sampled k points in Figs. 3(b) (c) is 8 with 
the Monkhorst-Pack method for the integration over the 
Brillouin zone^^, where the bands are fitted to sinusoidal 
forms and the tetrahedron method is employed. The re- 
sult in fact changes little when the number is increased 
to 18. 

We can see in the energy-resolved LDOS (Fig.3(b)) 
that, while the LDOS at Ef is small, LDOS is large 
around E Ep — ^ eV, which originates from 3(i levels of 
Ag. If we turn to the LDOS at Ep in Fig. 3(c), a notable 



peak is seen at the outermost oxygen, which suggests 
that the surface is in fact metalhc. This is a tail (in real 
space) of the bands originating from the MgO (marked 
in k space in Fig. 3 (a)). We have obtained a similar result 
for the case of two MgO layers (not shown). This is quite 
consistent with the experimental result^^ that unrecon- 
structed polar surface becomes metallic if the MgO film 
is ^ 5 MgO layers, and also with theoretical studies^"^. 
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FIG. 3. The band structure (a), the local density of states 
(LDOS) integrated over both x and y directions displayed as a 
function of energy and z{l. interface) (b), and the LDOS just 
at Ef as a function of z (c). The red mark in (a) indicates 
MgO- charactered bands, while vertical dotted lines in (c) rep- 
resent the optimized positions of atoms, with the separation 
between the outermost O and Mg highlighted by arrows. 
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FIG. 4. A plot similar to Fig.3 for one MgO(lll) layer. 



Now, for the thinnest possible (one) MgO layer in 
Fig. 4, we can see in the LDOS at Ef (c) that the peak at 
the O site continues to exist, which indicates that the one 
MgO (111) layer is still metallic as consistent with the ex- 
periment described above. To be more precise, the tail in 
real space still extends to the surface for the bands origi- 
nating from MgO although they [marked in Fig. 4(a)] are 
smaller in number for thinner films. The peak at the O 
site is seen to be much smaller in hight than in Fig. 3(c), 
so that the charge redistribution occurs less completely 
for one MgO layer than for three layers. Experimentally 



it is difficult to quantify the difference between one- and 
five-layer cases in the present EELS, but other techniques 
such as photoemission spectroscopy or energy loss spec- 
troscopy with higher energy resolutions should detect the 
difference. 

So a small electric field may remain for the incom- 
plete charge redistribution, but we will have to imple- 
ment a theoretical method to quantify the electric field^'^. 
Qualitatively, however, we can understand the incom- 
plete charge redistribution for one MgO layer as follows. 
We first note in Fig. 4(c) that the outermost O atoms sink 
deep into the Mg layer. Namely, the optimized separa- 
tion, in 2;-direction between the O and Mg atoms on 
the surface is 0.50 A, which is much shorter than that for 
the case of three MgO layers (with d = 0.90 A). This sug- 
gests that the electric field induced by the surface charge 
is small for one MgO layer. In other words, the energy 
loss due to the lattice distortion (Al) and that due to 
the buildup of the electric field (A^;) are relatively mod- 
erate for one MgO layer, so that the energy cost due to 
the charge redistribution (Ac) may exceed -\- Ae- 

This observation leads us to propose here an interest- 
ing phenomenon: if we can change the value of d which 
dominates Ai, + A^;, we can change the relative magni- 
tudes of Al -\- Ae and Ac- This implies that we should 
be able to induce a metal-insulator transition for the po- 
lar surface by controlling d. One of the simplest ways 
to change the value of d is to vary the lattice constant of 
the substrate. Namely, if we employ a substrate having a 
larger lattice constant, the separation between the hollow 
sites becomes larger, which in turn increases the distance 
between Mg atoms in x?/-plane. This will make O atoms 
sink more deeply into the Mg layer, and a smaller d will 
favor an insulating surface. 
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FIG. 5. A plot similar to Fig.3 for one MgO layer with a 
larger lattice constant (4.23 A) of Ag than that in Fig. 4. 

To confirm that the MgO surface becomes insulating 
for substrates having large lattice constants, we have cal- 
culated the band structure and LDOS for the interface 
system where we artificially make the lattice constant 
larger (4.23A Fig. 5) or smaller (3.90A Fig. 6). We can see 
that the larger lattice constant indeed makes the LDOS 



around Ep at the outermost O (Fig. 5) notably smaller 
than in Fig. 4, while the quantity becomes larger for the 
compressed substrate (Fig. 6). This shows that the sys- 
tem indeed resides in the vicinity of the metal-insulator 
transition, so it should be feasible to control the transi- 
tion. 
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FIG. 6. A plot similar to Fig. 3 for one MgO layer with a 
smaller lattice constant (3.90 A) of Ag than that in Fig. 4. 



From the behavior of the LDOS at Ep we can look 
into an even more fascinating phenomenon, i.e., the pos- 
sibility of a magnetic instability for the polar surface. As 
we have shown, the LDOS at Ef is large for the outer- 
most oxygen layer in thicker MgO films, where O atoms 
do not sink into the Mg layer. When the LDOS is large 
enough, we may expect a ferromagnetic instability. In 
fact, Goniakowskiet al.^ have already shown, by means 
of a spin density functional calculation for a five atomic- 
layer slab (0/Mg/O/Mg/O with unbalanced numbers of 
cations and anions), that the unreconstructed polar sur- 
face of MgO has a ferromagnetic instability. However, 
the structure of MgO (especially the value of d) is not 
optimized in that study, so the LDOS at Ej? should be 
overestimated. 
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FIG. 7. A plot similar to Fig. 3, where the spin density func- 
tional theory with GGA is adopted to calculate spin-resolved 
band structures. Blue (red) lines in (a,c) represent the ma- 
jority (minority) spin. 



So we have performed a spin density functional cal- 



culation for one-layer and three-layer MgO (111) on Ag 
with the optimization for d to study the magnetic insta- 
bility. The result shows that the three-layer system does 
have a ferromagnetic ground state, although the magne- 
tization (the difference between the numbers of opposite 
spins per unit cell) is small (~ 0.1). Figure 7 displays the 
band structure and the LDOS, where an exchange split- 
ting (~ 0.5 eV) is seen for the bands lying ~ 2 eV below 
Ep. One MgO (111) layer is paramagnetic as expected. 

One thing we have to note is that the unit cell in the 
present work as well as in Ref.^ is 1x1 in the x?/-plane, 
so that it cannot describe other magnetic states such as 
antiferromagnetic ones. Thus a calculation with larger 
unit cells will be required to really confirm whether the 
ground state can become magnetic, which is an impor- 
tant future problem. 
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